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ABSTRACT 

Although there has yet been no undisputed discovery of a stih-forming planet embedded in a gaseous 
protoplanetary disk, the cleared inner holes of transitional disks may be signposts of young planets. 
Here we show that the subset of accreting transitional disks with wide, optically thin inner holes 
of 15 AU or more can only be sculpted by multiple planets orbiting inside each hole. Multiplanet 
systems provide two key ingredients for explaining the origins of transitional disks. First, multiple 
planets can clear wide inner holes where single planets open only narrow gaps. Second, the confined, 
non-axisymmetric accretion flows produced by multiple planets provide a way for an arbitrary amount 
of mass transfer to occur through an apparently optically thin hole without over-producing infrared 
excess flux. Rather than assuming the gas and dust in the hole are evenly and axisymmetrically 
distributed, one can construct an inner hole with apparently optically thin infrared fluxes by covering 
a macroscopic fraction of the hole's surface area with locally optically thick tidal tails. We also 
establish that other clearing mechanisms, such as photoevaporation, cannot explain our subset of 
accreting transitional disks with wide holes. Transitional disks are therefore high-value targets for 
observational searches for young planetary systems. 

Subject headings: planet-disk interactions — protoplanetary disks — planets and satellites: formation 
— stars: pre-main-sequence — hydrodynamics 



1. INTRODUCTION 

Protoplanetary disks of gas and dust surrounding 
young stars are the birthplaces of planets — yet to date, 
there is no undisputed observation of a still-forming 
planet in a gaseous disk. Astronomers have there- 
fore paid considerable attention to transitional disks — 
protoplanetary disks with inner clearings — because such 
clearings occur as the result of tidal interactions between 
disks and planets (e.g. Lin & Papaloizou 1993). To 
date, over 30 transitional disks have been discovered with 
the IRS instrument aboard the Spitzer Space Telescope 
(MuzeroUe et al. 2010). If theorists were to demonstrate 
a planetary origin for any inner hole in a transitional disk, 
' that disk would harbor extraordinary discovery poten- 
tial for young, still-forming planets. Indeed, gaps in de- 
bris disks surrounding main-sequence stars harbor the gi- 
ant planets Fomalhaut b (K alas et al. 2008,), e Eridani b 
1 jHatzes ct al. 2000), the HR 8799 system (IMarois et al.l 
■ I2b0& ,2010, ) and /3 Pictoris b ()Lagrange et al.ll2010[ ). 

Inner holes in gaseous protoplanetary disks with ages 
^1 — 10 Myr are initially identified by the fiux deficits 
they leave in the near-IR and/or mid-IR part of the disk's 
spectral energy distribution (SED). The observed flux 
deflcits correspond to "missing" dust with characteristic 
temperatures of ~ 200 - 1000 K (e.glCalvet et al.| [2002l: 
IFurlan et al.l[200l IBrown et al.l 120071) . Model SEDs of 
the first transitional disks to be discov ered, GM Aur 
(iKoerner e t al.' 19931; iCalvet et all 12005") and TW Hya 
(jCalvet et a l. 2002). were based on disks with small in- 
ner holes only 3-4 AU wide. Soon after, investigations 
of planet-disk interactions demonstrated that such clear- 
ings co uld naturally be produced by a single J upiter-mass 
planet ()Ouillen et al.1 120041: iRice et al.ll2006D . Although 
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the planet opened only a narrow gap, not a hole ex- 
tending all the way to the central star, the dynamical 
studies showed that the residual gas and dust inside the 
gap would accrete onto the central star on the viscous 
timescale, widening the initial gap into a small, optically 
thin hole. 

The first hole-size measurements performed with mil- 
limeter interferometers s howed remarkable consistency 
with SED-based models (jBrown et al.ll2009() . However, 
further spatially resolved imaging of transitional disks 
brought a surprise. Some disks for which small, 3- 
4 AU inner holes were inferred based on SED modeling — 
including the canonical transitional disk, GM Aur — were 
shown to have large inner clearings of 15 AU or more 
once imaged at high angular re solution (jHughes et al.l 
2009HAndrews et al.ll2010l [20Tll). Y e t, as the models of 
Rice et al.l (|2006( ) and IQuillen et al.l (|2004[ ) hinted, each 
single planet can open only a narrow gap, which must 
then widen into a hole by viscous clearing of the ma- 
terial inside the gap. The extraordinarily long viscous 
timescales at distances of 15 AU or more meant the sim- 
ple theory that each transitional disk hole contained a 
single giant planet was no longer satisfactory. 

The prosaic explanation for wide holes in transitional 
disks — tidal truncation of the disk b y a stellar com- 
panio n, as observed in CoKu Tau/4 (jlreland fc KrausI 
I2008f ) — is also not likely responsible for the unusual and 
complex disks that concern us in this article. Based on 
observations with the Keck interferometer, IPott et all 
(|2010f) rule out stellar companions in 0.35-4 AU orbits 
at the 99.7% confidence level in the holes of DM Tau, 
GM Aur, LkCa 15, UX Tau A and RY Tau— at least 
three of wh ich are examples of t ra nsitional disks with 
wide holes (| Andrews et al.l 1201 ID . IKraus et al.) ()2011D 
also found no stellar companions in the inner holes of 
DM Tau, GM Aur, LkCa 15 and UX Tau A in their 
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non-redundant aperture masking survey of the Taurus 
star-forming region. The high accretion rates in these 
sources (see Table [1] and references therein) as well as 
the presence of inner disk dust and gas (e.g. iSalvk et"ari 
[2009':'EsD aillat et al.l[20Tnh . are also consistent with sub- 
stellar rather than stellar companions, as they indicate 
that disk truncation is incomplete. The set of transi- 
tional disks with large inner holes, up to 75 AU, and sig- 
nificant accretion represents a significant subset of tran- 
sitional disks, and a subset that continues to increase in 
number (Andrews ct al. 2011). In Table [U we list a sam- 
ple of eight rapidly accreting transitional disks that have 
holes between 15 AU and > 70 AU wide. 

In general, transitional disks have a variety of 
observational characteristics and may not repre- 
sent a single physical phenomenon. For example, 
clearing driven by the mag netorotational instability 
(jChiang fc Murrav-Clavl I2007f ) can widen dynamically 
cleared holes once such holes extend all the way to the 
dust-sublimation radius. In i )3.3l we explain in detail 
why magnetorotational clearing cannot be responsible 
for the particular set of transitional disks listed in Ta- 
ble [TJ Photoevaporation can also produce disks with in- 
ner cleared regions. Prior work on UV photoevaporation 
showed that evaporation-induced inner holes were incon- 
sistent with accretion rates hig her than ^ 10~^°Mq yr ~^ 
(| Alexander fc Armitagel[2007[ ). In essence, hole forma- 
tion can only occur once evaporation rates exceed accre- 
tion rates, as otherwise the hole will be replenished by 
accretion. More recently, however, it has been argued 
that X-ray photoevaporation should be the dominant 
disk-clearing proces s and produce rnuch h igher mass-loss 
rates than UV fe.g. lErcolano et al.|[2008[) . Thus, the set 
of parameter space consistent with photoevaporation has 
grown to include somewhat higher accretion rates. It is 
entirely possible that photoevaporation produces some 
transitional disks, while dynamical clearing by planets 
produces others. However, even X-ray photoevaporation 
models do not appear to be able to explain the existence 
of disks with both large holes > 15 -20 AU and signifi cant 
accretion rates > IQ-^-^Mq yr"^ (jOwen et al.ll2011| ). 

There is yet one more possible explanation for accret- 
ing transitional disks with wide holes — clearing by multi- 
ple planets, whose individual gaps could overlap to form 
a wide hole. Since planets have gap-crossing tidal tails 
and don't completely truncate the disk, a multiplanet 
system would allow gas and dust to accrete through the 
gap. In this work, we demonstrate that the particular set 
of transitional disks with wide holes and high accretion 
rates listed in Table [T] presents unambiguou s evidence 
for tidal clearing by multiple planets. While iZhu et all 
IJOll) found that wide, optically thin gaps in accreting 
disks could not be produced by planets alone, their con- 
clusions were based on the assumption of axisymmetric 
gas and dust in the hole. Instead, by taking into account 
the non-axisymmetric flow patterns produced by planets, 
we argue not only that multiple planets are consistent 
with high accretion rates and reduced infrared fluxes, 
but that they represent the only possible explanation for 
our sample of transitional disks. 

We begin our argument for planetary influence on tran- 
sitional disks in ^ where we demonstrate that plan- 
etary systems provide a viable explanation for the gas 
confinement, non-axisymmetry and mid-IR flux deficits 



observed in transitional disks with large holes. In ^ we 
show that the simplest planetary configuration — a single 
gas giant inside the hole — cannot explain the transitional 
disks with wide holes. We present dynamical models of 
multiple planets opening overlapping gaps in 21 In f}5]we 
use the density distributions predicted by the dynamical 
model in ^to produce SEDs, and discuss compatibility 
with observed photometry. In fj6] we delve into the pho- 
toevaporation mechanism in more detail, demonstrating 
conclusively that it is not a viable originator of our sam- 
ple of transitional disks. Finally, we present conclusions 
and testable predictions in |j7l 

2. PROPERTIES OF TRANSITIONAL DISKS WITH LARGE 

HOLES 

The primary observational characteristic that defines 
all transitional disks is that the near- and/or mid-IR 
fluxes are lower than those observed from 'typical' pro- 
toplanetary disks. Since the total flux emitted from an 
optically thick disk is not affected by small changes in 
the disk surface density, the deficit of flux in transitional 
disks is generally interpreted as requiring depletion suf- 
ficient to make the disk vertically optically thin at the 
flux-deficient wavelengths. However, it is important to 
keep in mind that the actual observable is the emitted 
flux, not the optical depth, t^. Knowledge of optical 
depth is required to convert observed fluxes into disk 
surface densities, but the assumption that is uniform 
throughout the hole, or even constant within a given an- 
nulus, can lead to serious errors in estimating the mass 
distribution in the inner disk. In particular, confined, 
optically thick 'streamers' within a hole can produce the 
same emergent flux as azimuthally symmetric, optically 
thin dust, while hiding a significant amount of the dust 
and gas that must move through the hole to match ob- 
served accretion rates. If the grains in a narrow tidal 
tail cover more surface area than the tidal tail itself, so 
that the tail is locally optically thick, some of the grain 
mass is effectively hidden. Clumping in Saturn's rings 
produces exactly this type of plateau in apparent optical 
depth as a function of particle surface density (Robbins 
et al. 2010). 

In the case of optically thick tidal tails surrounded by 
empty space, one can simply replace t^, with a geomet- 
rical "filling factor", /(r) — the fraction of an annulus's 
area covered by tidal tails — in the standard equation 
for disk flux: 



2TTrdrB^ [T (r)] f (r) / d"^ 



(1) 



where F^, is the flux at a given frequency, rjn and Tout 
are the disk inner and outer radii, respectively, is the 
Planck function in frequency space, d is the distance to 
the system, and T is the dust temperature at a given 
radius. An important difference between / and optical 
depth T^, however, is that / is purely geometric, and 
so independent of i^, while t,^ depends on the opacity k^, 
which is a function of v. Therefore, the mass distribution 
and/or temperature structure must be different in the 
two cases in order to produce the same spectrum. In 
addition, in a realistic disk, one must account for the 
optically thin disk atmosphere, which produces silicate 
emission features dependent on k^. 
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As an example of the ambiguities inherent in model- 
ing transitional disk SEDs, in Figure 1, we show two 
alternative inner disk dust distributions that produce 
similar near-IR spectra. To model the resultant flux of 
the models, we use d the radi ative transfer code RADMC 
(similar tO iDullemond fc Dominik 2004 . but without ver- 
tical structure integration) but optically thin compo- 
nents were also benchmarked with an analytical model. 
The first model consists of an optically thin region from 
0.2^1-1 AU (with Sdust = Eo(i?/AU)-°-S— used for this 
and all other models — and Eg = 5 x 10~* g cm^^) and 
a second optically thin region from 16-20 AU (with Sq = 
10"'' g cm~^). The second model consists of an optically 
thin region from 0.22-1 AU (with Eq = 3 x 10"^ g cm'^) 
and an optically thick streamer from 0.7 — 13 AU (with 
/ = 0.03 and H/R = 0.08 x {R/PMf/'^). Both models 
use the same standard Milky Way dust composition of 
80% small circums tellar silicates and 20% small carbon 
grains JWcingartn er fc Draini 120011: iLTfc Dr ainc 2001; 
iDraine et a l. 2007), with o paciti es from lOssenkopf et al] 
1199^ and Zubko eFall (|T996l) ("AGAR"). The car- 
bon grains, being featureless, only act to reduce the 
line/ continuum contrast, and so their quantity is degen- 
erate with other parameters, such as the inner disk small 
dust grain mass. Optically thin surface densities were 
chosen to roughly match GM Aur photometric data (al- 
though an exhaustive analysis of parameter space was 
not performed, and clearly these models have some trou- 
ble accurately fitting the observed silicate feature). The 
two models produce nearly identical emergent spectra. 
Figure 1 demonstrates that significant radial mass trans- 
fer can occur through transitional disk holes and still 
match photometric data, provided the mass is channeled 
into narrow streams such as planetary tidal tails. 

Besides extremely wide inner holes and high accretion 
rates, many transitional disks contain additional signa- 
tures of dyn amical perturbation by planetary systems. 
I Brown et al.l ( |2009j noted significant azimuthal asym- 
metries in the outer disk d ust distribution for LkHa 
330 and HD 135344 B, while IPontoppidan et al.l (pOOl 
noted asymmetries in the disks of HD 135344 B and SR 
21 on the scale of a few AU or lesfl Near-IR spectra 
were suggestive of radial dust asymmetries such as ec- 
centric inner holes in even early modeling studies (e.g. 
iGalvet et al.ll200"5l ). and recent studies show that both 
pre-transitional and transitional disks require multiple 
inner disk components, including nearly evacuated re- 
gions, in order to match observed flux e s and interfer- 
ometric visibilities ([Espaillat et all 120101 : [Andrews et al.l 
[20Tlt lAkeson et al .l l201lD. Similarly, SR 21 has a gas disk 
truncated at ~7 AU (jPontoppidan et al.l I2008D. while 
the dust disk is truncated at > 30 AU (jBrown et al.l 
120091: [Andrews et al.ll20TTI ). Finally, GO rovibrational 
lineshapes suggest that hot GO gas in transitional disks 
may be dynamically confined to smaller regions than in 
disks without inner holes (Salyk et al. 2011). The tidal 
barriers of massive planets would provide an extremely 
effective way to limit the radial distribution of hot gas in 

^ The measured inner radi us of GM Aur's dis k from near-IR 
interferometric measurements IIAkeson et al.|[2005l ). 

^ Although the asymmetries are observed in the astrometric pro- 
file of a CO emission line, the asymmetry can be present in either 
the gas or dust component. 



transitional disk holes. 

A final, unusual characteristic of transitional disks is 
that the y do not seern to ari se at any particular star age. 
Rather, iFurlan et al.l (|200 9l) show that the transitional 
disk fraction is 3-4% in the Taurus, Ghamaeleon and 
Ophiuchus star-forming regions, despite an age spread 
of 1.5 Myr betweeen the clusters. The transitional disk 
fractions reported by Gurrie & Sicilia-Aguilar (2011) also 
indicate that disks with inner holes (as opposed to ho- 
mologously depleted disks) are not preferentially found in 
clusters of a particular age. The inner-hole configuration, 
then, does not appear to represent an evolutionary stage 
all disks must go through — in which case there would be 
a pe ak in the tran s itiona l disk fraction for a given cluster 
age. iFurlan et al.l (|2009t ) use their observed transitional 
disk fractions to argue that the optically thin inner-hole 
configuration lasts about 100,000 years, which is simi- 
lar to Type II planet migr ation timescales from the 10- 
15 AU region (|Wardll 19971 ). Note, however, the warning 
by Gurrie & Sicilia-Aguilar (2011) that there is a degen- 
eracy between the time it takes a process to operate, and 
the frequency with which it operates. It is possible that 
10^ years is the hole clearing timescale rather than the 
lifetime of the inner hole. 

Planets, therefore, have the potential to provide all 
the unusual characteristics we observe in the subset of 
transitional disks under investigation here (Table [Ij : 

• Wide, apparently optically thin inner holes with 
r > 15 AU, 

• Near and mid-IR flux decrements as compared to 
classical protoplanetary disks, 

• High accretion rates of at least ~ 10~^ Mq yr~^ , 

• Gonfined, non-axisymmetrically distributed gas 
and dust. 

In the next section, we begin by examining the suitability 
of the simplest planetary configuration, a single gas giant, 
for opening the wide holes observed in our sample of 
transitional disks. 

3. ON THE FEASIBILITY OF A SINGLE-PLANET ORIGIN 
FOR LARGE HOLES 

Glearing by a single giant planet is a commonly invoked 
explanation for the gaps in transitional disks (e.g. Rice 
et al. 2003, Espafllat et al. 2008, Furlan et al. 2009). 
In ij3.1l we quantify the maximum gap width a planet 
can open, showing that it is not enough to explain our 
sample of accreting transitional disks with wide holes. In 
tj3.2l and 13.31 we show that even when one considers the 
combined effects of a gap-opening planet and accretion or 
spiral wake-driven clearing of material inside the gap, the 
timescale for opening a wide hole is outside the inferred 
10^-year transitional disk lifetime/clearing time. 

3.1. Maximum gap width 

To build a simple analytical framework for understand- 
ing the scope of planetary influence on disks, we begin by 
imagining the disk as a swarm of massless test particles. 
Treating the disk as if it were composed of test parti- 
cles on ballistic trajectories gives the maximum possible 
gap width a planet can open because gas pressure and 
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viscosity act to fill in the dynamically cleared gap. We 
estimate the maximum gap width by examining the time 
it takes for the planet to substantially alter the angular 
momentum of close-passing particles. The planet opens 
a gap when repeated encounters cause the particles at a 
given location to be strongly scattered within the age of 
the star/disk system. We define strong scattering as a 
change in the test particle's specific angular momentum 
of a factor of 2-3. The analytical framework used here 
applies to non-resonant particle orbits. 
We consider a test particle with Keplerian frequency 
and unperturbed orbital radius during a close en- 
counter with the planet, which has Keplerian frequency 
w and orbital radius rp. Both the planet and the test 
particle are initially on circular orbits. In a single en- 
counter, the particle's expected change in specific angu- 
lar momentum (Lin & Papaloizou 1993) is 



64 G^M^ro 
243 ^3^5 



2Kn 



1 2 



(2) 



and 



where Mp is the mass of the planet, A — \ro 
Kq and Ki are modified Bessel functions (Goldreich & 
Tremaine 1982). The rate of change in the particle's 
specific angular momentum due to repeated encounters 
with the planet is therefore 



■ _ She 
St ' 

where the time between encounters is 

St^ 2- 



(3) 



(4) 



The timescale over which the particle's specific angular 
momentum changes by of order unity is simply 



_h _rln 

^h^ h ' 



(5) 



Figure [2] shows the scattering timescale of disk par- 
ticles orbiting near two hypothetical planets, one with 
Jupiter's mass and the other with Saturn's mass. Both 
planets are placed 10 AU from a Sun-like star. Each 
planet's tidal radius, or Hill radius, is defined by 



Rh = Tp 



Mp 



1/3 



(6) 



given a stellar mass of A/* . For test particles within three 
Hill radii of the planet's orbit, the scattering timescale 
is 10^ years or less. Tidal interactions with planets are 
therefore extremely efficient at opening gaps in proto- 
planetary disks. Notably, Crida et al. (2006) found that 
even in strongly accreting disks with kinematic viscosity 
of order 10^^ cm^ s~^ at 1 AU, corresponding to an accre- 
tion rate of « 5 x 10^^ Af© yr"^ (Dodson-Robinson 
et al. 2009a), a Jupiter-mass planet (IMj) can open a 
gap in only 1000 orbits. Encouragingly, planets can open 
gaps even in disks where strong accretion would readily 
fill in gaps created by photoevaporative winds. 

There is, however, a limited scope to the tidal in- 
fluence of the planet. As we see from Figure [U the 
scattering timescale for particles located four Hill radii 
or more from the planet approaches the ages of transi- 
tional disk host stars. Setting a high scattering timescale 



limit of 5 X 10^ years, which is one-half to one-quarter 
the esti mated stellar age in the Taurus- Auriga asso- 
ciation (|Palla fc Stahleii |2000[ ) and five times higher 
than the inferred tr ansitional disk lifetime/ clearing time 
()Furlan et al.l[200l . we find that the gap opened by a 
planet may not exceed ^ 5 Hill radii. And, since the Hill 
radius scales as Mp^'^, simply increasing the planet mass 
does not provide a commensurate widening of the gap. 
Even for a 5Mj planet at 10 AU, the Hill radius is only 
1.2 AU, which allows a gap to open from 4 AU-16 AU. 
This gap is simply not wide enough to account for the 
observed hole sizes of any of the transitional disks listed 
in Table [1] 

Given the sharp increase in scattering timescale as a 
function of distance from the planet (Figure [2]), our max- 
imum gap width is not highly sensitive to our assumed 
gap-opening timescale. Our order-of-magnitude estimate 
of maximum gap width agrees with the numerical work of 
Crida et al. (2006). Using analytical and numerical cal- 
culations of the balance between torques from the planet 
and viscous stresses, they found that a Jupiter-mass 
planet in a disk with aspect ratio H/r = 0.05 opened 
a gap with a full radius of 5.5 Rh between the planet 
orbit and the primordial, unperturbed disk. The width 
of the optically thin gap is, of course, much smaller than 
the full center-to- wall width. Our calculation of 5 Rh as 
the maximum gap size is actually rather generous. 

3.2. Interior disk clearing by viscous accretion 

Our calculations demonstrate that a single planet can- 
not clear a gap wider than ^ 5 tidal radii over typical 
T-Tauri star ages. A Jupiter-mass planet at 10 AU could 
therefore open a gap extending from 6.5-13.5 AU. How- 
ever, the star is accreting: could the accretion flow drain 
the disk interior to the planet's orbit (which we hereby 
label the "interior disk" ) within the transitional disk life- 
time and/or clearing time? For a disk with steady-state 
mass transfer, the diffusion timescale is 



^diff 



(7) 



where Af* is the star accretion rate and M{rg) is the 
interior disk mass. A derivation of Equation [7] is in Ap- 
pendixEl For narrow annuli of material located near the 
star, so that < 2 AU, clearing of the interior disk by 
viscous accretion seems possible. Consider a disk with a 
surface density profile of 



E(r) = (808 g cm" 



1 AU 



(8) 



and an inner edge due to magnetospheric truncation at 
0.1 AU. This is the pre-clearing density profile of the 
CM Aur disk constructed by Hughes et al. (2009), who 
inferred a total disk mass of O.lGAf© and an outer radius 
of 300 AU from millimeter imaging. If the star accretes at 
a rate of 10~*A/q yr^^, the diffusion timescale for mate- 
rial in the inner 2 AU is 10^ years, similar to transitional 
disk lifetimes/clearing times (Furlan et al. 2009). 

However, what about holes that are 20 AU or larger? 
Once can place a planet at most bRn from the desired 
edge of the hole, or at about 14 AU for a 3Mj planet or- 
biting a Sun-like star. (Since Rh scales linearly with the 
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planet's orbital radius, we open the widest gap if we place 
the planet as far from the star as possible.) The inside 
edge of the planet's gap is then at about 8 AU. For the 
inferred GM Aur density profile, the diffusion timescale 
for material at 8 AU is now 4x10^ years, edging toward a 
significant fraction of star ages in the Taurus, Ophiuchus 
and Coronet associations ()Furlan et al.l 120091 ) . Further- 
more, the diffusion timescale is a really an e- folding time 
(Appendix [X]) , not a time in which complete clearing of 
interior disk material occurs. 

To check how many e-foldings it takes to make our inte- 
rior disk optically thin at millimeter wavelengths, where 
the wide transitional disk holes have been imaged (Brown 
et al. 2008, Hughes et al. 2009, Andrews et al. 2011), we 
use the opacity law 



KA = 0.02(^1 cmV , 



(9) 



(jBeckwith et al.l[l990[ ). In Equation [9l k\ is the total 
mass opacity (gas -|- dust) at a given wavelength A, /? 
is an exponent near unity, and gas and micron and/or 
submicron grains are mixed with a mass ratio of 100:1. 
Throughout this article we will use the convention that 
K\ refers to an opacity based on total mass, gas and dust 
together, and n^dust gives an opacity referenced to the 
mass of dust grains only. 
The vertical optical depth of the disk is simply 



(10) 



A surface density of only 20 g cm^^ of gas and dust 
mixed with a 100:1 mass ratio is sufficient to make a 
disk optically thick, r = 1, at A = 0.4 mm. Diffusing 
the 100 g cm-2 initially at 8 AU in our reconstructed 
GM Aur disk down to optically thin levels would take 
two e-foldings, or 8 x 10^ yr. Even for our hypothetical 
gap with an inner edge at 2 AU, for which the charac- 
teristic diffusion timescale is seemingly reasonable, three 
e-folding times are required to get to optically thin levels, 
or 3 X 10^ yr. If a planet is far enough from its host star 
for the outer edge of its gap to be at 15-20 AU or more, 
there is a severe timescale problem with clearing its in- 
terior disk to optically thin levels by viscous accretion. 
As we explore in the next section, the mid-IR limits on 
leftover surface density from the interior disk are even 
stronger than the limits from millimeter observations. 

3.3. Interior disk clearing by other mechanisms 

Interior disk material could be cleared via the scenario 
proposed by Varniere et al. (2006), in which spiral waves 
excited by a single planet drive enhanced angular mo- 
mentum transport that clears the interior disk within 
^ 0.16idiff- Using 2-d hydrodynamic simulations of 
planet-disk interactions, Varniere et al. show that within 
2,000 orbits, the planet's spiral wakes reduce the interior 
disk surface density by about a factor of 10. To test the 
viability of spiral wake-driven clearing, we once again set 
a target hole size of 20 AU to match the GM Aur disk 
and place a hypothetical planet at 14 AU, 5Rh from 
the edge of the hole. Our hypothetical planet takes 10^ 
years to complete 2,000 orbits and open its incomplete 
hole — so far the clearing time matches the transitional 
disk lifetime predicted by Furlan et al. (2009). 



But an axisymmetric factor of 10 surface density reduc- 
tion is not nearly enough to make an optically thin hole in 
the mid-infrared. Using a micron/submicron grain opac- 
ity of KA.diist = 500 cm^ g~^ at 13 fxm (jOssenkopf et al.l 
Il992f ). a gas/small grain mass ratio of 100 and Equa- 
tion [TOl we see that the total (gas-f dust) surface den- 
sity required for a disk to be vertically optically thin at 
13 /im — the spectral index that usually diagnoses flux 
deficits in transitional disk SEDs (Furlan et al. 2006) — is 
only 0.2 g cm~^. For spiral waves to reduce the surface 
density from 100 g cm~^ to 10 g cm~^ at our 8 AU in- 
terior disk edge still leaves 50 times too much material 
to have an optically thin disk at 13 /im. Single plan- 
ets and their spiral wakes therefore do not clear enough 
material from the interior disk to explain observed transi- 
tional disk SEDs. However, our submicron grain opacity 
value is calculated assuming the grains have maximum 
sizes of 1 /im or less. It is possible that the combination 
of a single planet, wake-driven hole clearing, and lock- 
ing away 98% of the solid mass in pebbles, planetesimals 
and/or planet cores could be consistent with observed 
transitional disk fluxes. 

We have also considered whether photoevaporation 
or the magnetorotational instability ( "MRI" , Chiang & 
Murray-Clay 2007) could clear the interior disk material. 
As we discuss in detail in fj6l phoevaporative clearing is 
only efficient at radii > 3 AU. In such a scenario, the in- 
terior disk still clears via viscous accretion, and is subject 
to the time constrain ts we have just discussed . The MRI 
scenario proposed bv lChiang fc Murray-Clavl (|2007f ) also 
cannot clear the interior disk material. If the inner disk 
is already cleared of most small dust grains. X-rays can 
ionize material at the edge of the hole and produce inside- 
out hole widening. The accretion rate is then found to 
scale roughly as the square of the hole size. Inside-out 
MRI clearing is therefore not effective at initiating the 
clearing of material from the interior disk (ie, at small 
disk radii). 

Transitional disks such as GM Aur and LkHa 330, 
with inner hole radii of 20 AU and 50 AU, respectively 
(Hughes et al. 2009, Brown et al. 2007), cannot be ex- 
plained by a single orbiting planet. Even if one takes into 
account accretion or spiral wave-induced clearing of the 
material left inside the planet's gap, the inner hole size 
that can be cleared by a single planet in of order 10^ years 
is limited to a few AU. However, multiple planets or- 
biting the same star would open adjacent gaps which, 
combined, could form a large inner hole. In the next sec- 
tion, we use numerical simulations to demonstrate how 
the "overlapping gap effect" can explain even wide inner 
holes in accreting disks. 

4. DYNAMICAL CLEARING BY MULTIPLE PLANETS 

Our working picture of a transitional disk is an op- 
tically thick outer disk surrounding a wide inner hole 
that produces little near- and mid-infrared flux. The star 
must also be able to accrete material through the hole 
at a high rate, as much as 10~^Mq yr~^. Even though 
a single planet can only open a narrow gap, a system of 
multiple planets could produce a wide hole if the planet 
orbits were close enough together to create overlapping 
gaps. Because the GM Aur system is one of the best- 
characterized transitional disks known, we use its overall 
mass, hole size and accretion rate (Table [TJ as guidelines 
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for our study of multiple planets in transitional disks. 

We describe our hydrodynamic simulations of hole 
opening by interacting planets in ^A.ll Then, in H4.21 we 
show that the non-axisymmetric flow patterns produced 
by planetary systems are the best possible explanation 
for the high accretion rates of transitional-disk host stars. 

4.1. Hydrodynamic simulations 

Our planet-disk simulations were conducted using the 
publicly available version of FARGO, a polar, grid-based 
2-d hydrodynamic code designed specifically for planet- 
disk interactions written by Frederic Masset (Masset 
2000). Like all polar hydrodynamic codes built on the 
operator-splitting technique, FARGO uses a sequence of 
three steps — source step, radial transport step, and az- 
imuthal transport step — to evaluate the transport of any 
hydrodynamic quantity. The FARGO transport algo- 
rithm differs from standard codes in that the azimuthal 
advection is carried out in a reference frame that rotates 
at the local Keplerian speed, so that hydrodynamic quan- 
tities are advected using the residual velocity, — v^, 
instead of the inertial azimuthal velocity. A subsequent 
step corrects for the azimuthal shifts in the field in- 
troduced by subtracting out the Keplerian speed. The 
advantage of the FARGO transport algorithm is that 
timesteps are not Courant-limited by the short dynam- 
ical time at the inner disk boundary. FARGO produces 
up to a tenfold speedup in 2-d hydrodynamic calcula- 
tions over a traditional advection procedure, making sim- 
ulations that cover a significant fraction of the star age 
tractable. 

In order to test whether a stable multiplanet system 
could reside inside a transitional disk hole, our simula- 
tions include both migration-causing torques from the 
disk and dynamical interactions between the planets. Us- 
ing the 20 AU hole size of GM Aur's disk and our single- 
planet gap size of 5Rh as a guide, we set up a disk with 
three planets located about six mutual Hill radii apart. 
Chatterjee et al. (2008) defined the mutual Hill radius in 
planet encounters as 



V 3M, 



(11) 



Tightly-packed systems with planet orbits spaced less 
than ARHm apart strongly scatter on 1-10 Myr 
timescales and are not likely to survive intact (Chatter- 
jee et al. 2008, Scharf & Menou 2009). For planet orbits 
QRiim apart, the outer ARh boundary of the innermost 
planet, measured in individual Hill radii, coincides with 
the inner ARh boundary of the next planet in the system. 
Our orbital spacing therefore favors the creation of over- 
lapping gaps, which can reach 5Rh, yet maximizes the 
likelihood of long-term dynamical stability. To recreate 
the 20 AU hole in the GM Aur disk, we place three plan- 
ets, each of three Jupiter masses, at distances of 14.3 AU, 
6.3 AU and 2.7 AU. Our planetary configuration would 
also provide a reasonable match for the holes of DM Tau 
or DoAr 44 (19 AU and 30 AU, respectively; Andrews et 
al. 2011). 

The disk has a constant aspect ratio of H/r = 0.05, 
where H is the pressure scale height. We place the outer 
boundary at 38 AU and the inner boundary at 0.72 AU. 
Our evenly spaced grid contains 220 radial zones and 384 



azimuthal zones. The inner boundary is open so that 
mass can leave the grid, simulating accretion onto the 
star. We also allow mass to enter the grid at the outer 
edge of the disk, mimicking a steady-state mass transfer 
profile. The initial surface density profile of the disk is 
the same as in Equation [8] (fj3]). If we extrapolate the 
surface density profile out to 300 AU, the observed size 
of the GM Aur disk (Hughes et al. 2009), the overall disk 
mass would be 0.16Mq. 

Accretion through the disk is parameterized using the 
a- viscosity, where the turbulent diffusion coefficient is 



ac^H. 



(12) 



In Equation I12| Cs is the sound speed and a is the effi- 
ciency with which gravitational energy is extracted from 
the shear ffow and transformed into turbulence. Our 
choice of a was motivated by studies of turbulence caused 
by the magnetorotational instability or MRI, the leading 
contender for the viscosity source in protostellar disks 
(Balbus & Hawley 1991, Balbus et al. 1996). Turner & 
Sano (2008) find that a < 10~^ inside the optically thick, 
nonionized "dead zone" within two pressure scale heights 
of the disk midplane. We adopt a turbulent efficiency of 
a = 0.002, consistent with mass-weighted estimates of 
local accretion rates with turbulence driven by the MRI 
(Turner & Sano 2008, Turner & Drake 2009). All param- 
eters in the FARGO simulation are summarized in Table 

m 

Like Zhu et al. (2011), we allow FARGO to apply only 
the Type II migration torques. In Type I migration, 
which occurs before the planet opens a gap, the disk's 
outward pressure gradient amplifies the torques from the 
planet's wake at the outer Lindblad resonances, forcing 
the planet to move inward at an extremely high rate 
(Ward 1997, Tanaka et al. 2002). In the limit of low 
thermal diffusivity, the nonlinear torque caused by the 
"horseshoe" of gas and dust corotating with the planet 
saturates and the only net torque on the planet comes 
from the Lindblad resonances (Ward 1991). Yet a new 
series of investigations of planet migration in viscous 
disks with nonzero thermal diffusivity demonstrates that 
the corotation torque may halt inward migration or even 
drive outward migration (Paardekooper & Mellema 2006, 
Kley & Crida 2008, Kley et al. 2009, Paardekooper et al. 
2010, Paardekooper et al. 2011). Since FARGO does not 
treat heat or momentum diffusion, the migration rates 
calculated in the Type I limit — either no gap or an in- 
complete gap — may not accurately describe the motion 
of a planet in a physically realistic disk. To prevent im- 
proper computation of Type I migration torques, we first 
follow the gap-opening process for 200 orbits of the out- 
ermost planet (12,000 years) with the migration feature 
of the FARGO code turned off. We then "release" the 
planets, turning on the migration torques and the N- 
body integrator to track planet-planet interactions. 

Unlike Zhu et al. (2011), we do not include accretion 
onto the planets, which can be parameterized in FARGO 
by specifying the percentage of mass entering the planet's 
Hill radius that the planet can accrete, p. Over the 
course of both the initial 12,000-year gap-opening epoch 
in our simulations and the > 10^-year length of the hole- 
opening phase (Furlan et al. 2009), a young gas giant 
shrinks by at least two orders of magnitude, dramatically 
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decreasing its accretion efficiency. A physically realistic 
estimate of p would therefore require both an adaptive 
mesh to model the circumplanetary disk in detail and ro- 
bust knowledge of the interior structure and evolutionary 
stage of each planet. Because the 3Mj planets we have 
placed in our model transitional disk should be much 
smaller than their Hill spheres for most of the simula- 
tion, we have assumed that accretion onto the planets 
clears a negligible amount of disk material. In Appendix 
[B]we review the evolution of a young gas giant's interior 
structure in more detail and discuss the effects of these 
changes on the planet-growth process. 

4.2. Hole structure and appearance 

Figure |3] shows the surface density of gas and dust as a 
function of time in our simulation. A wide hole created 
by three overlapping individual gaps opens within about 
300 orbits of the outermost planet, or 1.8 x 10^ years. The 
color scale in the topmost row of Figure [3] is designed to 
mimic sub-mm opacities. Based on the opacity law in 
Equation |9l disk gas and dust mixed with a 100:1 mass 
ratio with surface density 20 g cm~^ is optically thick 
(r > 1) at 400 /im. All parts of the disk rendered in white 
are optically thin in the sub-mm. After 400 orbits of the 
outermost planet, or 2.4 x 10** years, the optically thin 
sub-mm hole extends to about 17 AU. At A = 1 mm, the 
hole extends to about 20 AU. The appearance of a wide 
hole consistent with sub-mm observations is not merely 
a result of having assumed a low-density disk: our ini- 
tial surface densities are entirely based on the mass of 
the GM Aur disk inferred from millimeter observations 
(Hughes et al. 2009), which is quite high at O.I6M0. 
Interestingly, Salyk et al. (2009) note that transitional 
disks with measurable accretion tend to have higher in- 
ferred masses than "classical" disks without holes, which 
supports planet formation as a likely explanation for ob- 
served hole structures (although the low-mass transition 
disk fraction is not yet well known; Andrews et al. 2011). 

Furthermore, all transitional disks have mid-IR flux 
deficits indicating their inner holes are optically thin 
from 10-25 /im. On the bottom row of Figure [31 we 
stretch the color scale to mimic 13 /im optical depth. At 
13 /im, micron and submi cron grains produce an opacity 
of i^\,dust ~ 500 cm^ g ^ ([Ossenkopf et al]|1992f ) that is 
not very sensitive to size distribution as long as the max- 
imum grain size is under 10 /im. The r = 1 threshold 
is then at a dust-only surface density of 0.002 g cm~^ 
of micron and submicron grains. Assuming a gas/solid 
mass ratio of 100 and making the assumption that 100% 
of the disk's solids are in micron and submicron grains 
(as opposed to pebbles, planetesimals or planet cores), 
the gas and dust together must have a surface density 
S < 0.2 g cm^^ to be vertically optically thin at 13 /im. 
Again, all optically thin regions of the disk in the small- 
grain approximation are rendered in white. Globally, the 
inner hole has a filling factor of / « 0.5 and is overall 
optically thin. 

Where a significant portion of the grains have formed 
pebbles, planetesimals or planets, the optically thin 
threshold surface density of gas and solids mixed with 
a 100:1 mass ratio is higher than 0.2 g cm~^. In Fig- 
ure m we once again render the 400-/im and 13-/im op- 
tical depths calculated from orbits 197 and 395 of our 
hydrodynamical simulation, accounting for the fact that 



a significant fraction of the available solid mass must 
be locked away in giant planet cores, residual planetesi- 
mals, or planetesimal precursor particles. Our fiducial 
opacity of Kx,dust = 500 cm^ g~^ for dust at 13 /im 
is valid to about 20% for any single grain size or size 
distribution as long as the maximum grain size present 
is of order 1 /zm. Ye t Mie scattering calculations by 
lOssenkopf et al.l ()1992D show that for the MRN gram size 
distribution, 

= -3.5, (13) 

a log a 

wher e n is the number o f grains and a is the grain ra- 
dius (|Mathis et al.l[T977h . significant 13-/im opacity re- 
ductions begin when the maximum grain size is extended 
to 10 /im. Extending the MRN distribution to to mil- 
limeter, "pebble" sizes produces a factor-of-10 opacity re- 
duction over our fiducia,! value of kx dust- 

In Figure 13] we 

assume that 80% of solids have either grown to millime- 
ter sizes or become incorporated in planetesimals/planets 
and the remaining 20% of the solids are micron and sub- 
micron grains which provide the observed fiux. The op- 
tically thin threshold surface density of gas and solids 
mixed with a 100:1 mass ratio is then E < 1 g cm~^ at 
13 /im and E < 100 g cm"^ at 400 /im. At 13 /im, al- 
most 90% of the inner hole's surface area contains only 
optically thin dust and the optically thick tidal streams 
cover only 10% of the hole, / w 0.1. 

In Figures [3] and m we see the detailed structure of each 
planet's tidal tail, as well as small pileups of material be- 
tween the planets. Far from being problematic sources 
of excess flux, these pileups and tails provide a critically 
important feature of our model: the ability to transport 
mass through an apparently optically thin hole. The 
apparent contradiction between optically thin holes and 
high stellar accretion rates comes from the assumption 
that gas and grains are evenly and axisymmetrically dis- 
tributed within the hole (Espaillat et al. 2008, Espaillat 
et al. 2010, Zhu et al. 2011). For uniformly distributed 
grains and a gas/solid mass ratio of 100, our 0.002 g cm^^ 
of dust, T ~ 1 threshold would allow less than eight lu- 
nar masses of dust in the hole, or O.IAIq. If 80% of the 
solid mass is in pebbles, planetesimals or planet cores 
and only the remaining 20% provides 13 /im opacity, the 
total mass of solids inside the hole could still not exceed 
0.5Afe. 

By contrast, after 24,000 years of evolution, our simu- 
lated disk's inner hole contains about O.95M0 of solids — 
yet this material covers between 10% and 50% of the 
surface area of the hole, depending on the assumed 
gas/small grain mass ratio. In our fargo simulation, 
the inner hole contains about Saturn's mass {95M,^) 
of gas and dust together. We have therefore used the 
non-uniform, non-axisymmetric property of flow patterns 
modulated by giant planets to at least double the mass 
of small grains consistent with an apparently optically 
thin inner hole. 

Another new piece of observational evidence supports 
the idea that gas in the inner holes of transitional disks 
is conflned by planets. Using rotation diagrams of the 
CO fundamental rovibrational ladder, Salyk et al. (2011) 
found that simple emission models are consistent with 
systematically lower surface areas for transitional disks 
than for gap-free protoplanetary disks. The tidal barriers 
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of planets in the inner holes of transitional disks would 
confine hot gas to tidal tails or narrow annuli between 
planets and therefore lead to a lower emitting area for 
hot gas in transitional disks. 

One issue highlighted by Zhu et al. (2011) was their dif- 
ficulty keeping the planetary system stable. They found 
that migration drives the planets into 2:1 resonances, 
which leads to strong scattering. The high migration 
rate in their simulations is due to their choice of a in the 
standard a-viscosity parameterization: Zhu et al. used 
a = 0.01, where we use a = 0.002, motivated by sim- 
ulations of MHD turbulence f §4.ip . Type II migration 
operates on the disk's viscous timescale, which is shorter 
for higher values of a. The planets in our simulations do 
begin to migrate after they are released, particularly the 
outermost planet, which moves from 14.3 AU to 10 AU 
over the first 5 x 10** years of the simulation. Yet our 
orbit integration shows that planetary system remains 
stable for more than 10^ years. Although physical under- 
standing of turbulent viscosity generation in disks is still 
limited, it is possible to construct a physically realistic 
transitional disk in which migration does not destabilize 
the planetary system. 

5. MODEL TRANSITIONAL DISK SEDS 

In this section we consider whether the multiplanet 
dynamical model shown in Figure [3] is consistent with 
observed photometry of transitional disks. Our goal is 
to discuss the degree of compatibility between our model 
and observed fluxes; an exact fit of SED model to data 
is beyond the scope of this work, and would be ham- 
pered by significant limitations. Among the limitations 
is the fact the dynamical model predicts only the gas 
distribution, while most observables are based on emis- 
sion from small dust grains. Because such grains would 
be affected by pressure gradients set up by the complex 
gas structure, and because the expected gas/small-dust- 
grain ratio is essentially unknown, the translation from 
gas distribution to an emission spectrum is not trivial. 

Nevertheless, we make some simple assumptions here 
and produce models designed to give us a qualitative 
understanding of the compatibility between the multi- 
planet simulation and transitional disk observations. In 
particular, we believe that the most stringent require- 
ment is that mid-IR fluxes be as low as observed, as most 
other features of the disk SED can be flt by tweaking any 
of the numerous (and degenerate) parameters that enter 
sophisticated disk models. We begin by comparing our 
model to photometric data of the GM Aur system. How- 
ever, we find that our model is in some ways more natu- 
rally consistent with pre-transitional disks. We therefore 
also discuss the pre-transitional disk DoAr 44. 

Our simple model consists of three components: an 
optically thin inner disk from i?in to 0.7 AU, an op- 
tically thick component from 0.7-13 AU with a con- 
stant filling factor / meant to represent the dynamical 
model described in ^ and an optically thick outer disk 
from 13 AU to i?out (see Figure [5]). Model SEDs are 
created using the Monte Carlo radiati ve transfer code 
RADMC (Dul lemond fc Dominil] l2004f ). again without 
vertical structure integration. Stellar pa rameters and 
disk i nclination are set to known values (jHughes et al.l 
[200l [Andrews et aLllMl . We use a standard Milky 
Way dust composition of 20% carbon and 80% small 



(sub- j itm) circumstellar silicates (iWeingartner fc Draind 
[200ltlLi fc Drainell2001l:rDraine et al.ll2007D . i?in was set 
at the measured dust inner radius, if known, or at roughly 
the dust sublimation radius, and the mass and radial 
extent of optically thin dust were adjusted to roughly 
match observations. To simulate an outer disk that is 
not directly irradiated by stellar flux but instead mostly 
shadowed by tidal streams, we took the difference be- 
tween a model from it!in to -Rout and another model from 
i?in to 13 AU — although, as is apparent in Figure EJ the 
observations do seem to require p artial illumination to 
'puff up' the hole's outer edge (e.g. [Andrews et alll2011l : 
iHughes et al.l 120091) . However, the flux from the outer 
disk is not the focus of this work, and we have made no 
effort to fully explore parameter space for these compo- 
nents. 

Instead, we are concerned with the near- and mid-IR 
flux contribution from the confined tidal streams in our 
dynamical model. We begin by discussing the mid-IR 
fluxes, which originate from the ~ few AU region of 
the disk. The left-hand side of Figure [5] shows the az- 
imuthally averaged dust surface density for our dynami- 
cal model after 400 orbits of the outermost planet. If one 
assumes the mass is spread evenly over all azimuths the 
surface density is almost everywhere optically thick at 13 
(assuming an opacity of Kx,dust = 500 cm^ g~^ for 
dust grains). Not surprisingly, since transitional disks 
by definition have very low 13 /im fluxes, an optically 
thick disk model with / = 1 (red squares) is not consis- 
tent with observations. In reality, however, the density 
distribution predicted by the dynamical model is not ev- 
erywhere optically thick, as the mass is shepherded due 
to the presence of the planets. In the language of fJS] the 
azimuthal flUing factor, /(r), (gray solid line), ranges 
from ^0.1-1 if the gas/solid mass ratio is 100 and all 
solids are in micron and submicron grains. Extending the 
grain-size distribution up to millimeter radii or locking 
away solids in planetesimals and planets correspondingly 
lowers / (see Figures [3] and |4]) . 

Full 2-d modeling of the emergent flux from the dis- 
tribution predicted by our hydrodynamical simulation is 
beyond the scope of this analysis; instead, we make the 
simple assumption that the flux scales proportionally to 
some characteristic average value of /, and ask whether 
the resulting fluxes are consistent with photometric data. 
We implicitly assume that optically thin parts of the disk 
are completely devoid of grains, and so this can be con- 
sidered a 'best-case' scenario. In reality, the optically 
thin regions would need to have an average r lower than 
that derived from axisymmetric models. On the other 
hand, we begin with a conservative assumption about 
the gas/small-dust- grain ratio, and later discuss the im- 
plications of this choice. 

The middle panel of Figure [5] compares model SEDs 
with photometry and IRS spectroscopy of GM Aur. The 
blue curve (triangles) in Figure [5] shows the flux pro- 
duced if the inner hole has / — 0.5, a representative 
value for our dynamical model if the gas/dust mass ratio 
is 100. Although the flux is signiflcantly reduced from 
the optically thick case, fluxes are still inconsistent with 
observations. Instead, to match observations, / must be 
signiflca ntly smaller, o r ^0 .03 (black curve). For com- 
parison, iCalvet et all (|2005D estimate r ~ 0.002 at 10 
/xm. Note that we have not made an effort to fit the 
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silicate feature exactly; however, one can add additional 
contrast by adding small regions of optically thin dust in 
the inner disk, and/or by adjusting the grain composi- 
tion (and thus opacity). Pre-transitional disks like DoAr 
44 (right hand side of Figure [5]) are similarly inconsistent 
with the completely optically thick model (red squares), 
as well as the / = 0.5 case (blue triangles), using stan- 
dard assumptions. However, since the mid-IR fluxes in 
this case are significantly higher, the required reduction 
in flux is more moderate — close to / = 0.1 (black curve). 

Are the observations consistent with the dynamical 
model? As discussed in §4.2) / is essentially anti- 
proportional to the fraction of solid mass locked in peb- 
bles, planetesimals and planet cores, so one way to match 
observations of a disk like GM Aur is simply to reduce 
the mass fract ion of sma l l grain s in the inner disk, as also 
suggested by Zhu et al.l ()2011| ) . That much of the mass 
has been incorporated into large grains is certainly a rea- 
sonable possibility, given the fact that large planets have 
already formed. It is important to keep in mind that our 
simple SED models ignore the flux contribution from ar- 
eas where r < 1 locally, and so the reduction in flux with 
changes in / due to grain growth is optimistic. Nev- 
ertheless, with a power-law size distribution of silicate 
grains with a minimum grain size of amin = O.Q05/im, 
setting a„iax = 1 mm reduces 13 /xm opacity by a fac- 
tor of '^lO relative to a distribution composed entirely 
of small grains. Therefore, even moderate grain growth 
can easily make these models a reasonable possibility. 

As an alternative to flux reduction by grain growth, 
one could geometrically constrain the small grains, re- 
ducing their filling factor relative to that derived for the 
gas. It is well known that orbiting dust grains subject 
to a gas pressure gradient can migrate towards pressure 
maxima (Weidenschilling 1972), and so it is not unrea- 
sonable to expect some additional shepherding of dust 
grains relative to the gas. 

We now consider the near-IR fluxes. One notable fea- 
ture of our dynamical model (see Figure [S]) is the per- 
sistence of a small annulus of optically thick material 
between the innermost planet and the star. The dynam- 
ical influence of the innermost planet does not extend 
to these radii, or to radii smaller than 0.7 AU, which 
is the inner boundary of the disk in our dynamical sim- 
ulation. Thus, the very inner disk is not dynamically 
cleared, and, if present, contributes significant near-IR 
flux. As discussed in ^ it is possible to hide some opti- 
cally thick, but very geometrically thin, material in the 
GM Aur inner disk. However, near-IR fluxes are cer- 
tainly inconsistent with the large optically thick dust rim 
extending from 0.2 to 1.5 AU that the dynamical model 
leaves behind. Therefore, for disks like GM Aur, this in- 
ner disk material must be cleared by another mechanism, 
such as viscous dissipation. Our system of planets was 
specifically designed to have a viscous clearing timescale 
for inner disk material of ~ 10^ years or less, so this 
possibility is within reason. 

For so-called pre-transitional disks such as DoAr 44, 
however, near-IR fluxes actually require a significant 
amount of material close to the dust sublimation radius 
(Espaillat et al. 201 0; see also Figure [5l). This interior 
pileup, also found bv lRice et al.l ()2006[ ). is a natural con- 
sequence of multi-planet systems of the type described 
here. In addition, although pre-transitional disks are es- 



sentially defined by the presence of this inner disk mate- 
rial, even "traditional" transitional disk sou rces typically 
have measurab le excess dus t emission (e.g. lEisner et al.l 
l2006t .Espaillat et al.l l2010f) as well as emission from 
hot CO ga s ( Rcttig ct aL 200j IPontoppidan et al.ll2008l: 
ISalvk eta!.. 2009) originating at r < 1 AU. Thus, our 
model suggests that the range of observed near-IR fluxes 
represents a range of stages of dissipation of the material 
interior to the innermost planet. 

6. OTHER HOLE-OPENING MECHANISMS: 
PHOTOEVAPORATION 

In 52] we showed that systems of multiple planets pro- 
vide a robust explanation for the observed characteristics 
of transitional disks: high M, wide holes that appear op- 
tically thin, and confinement of hot gas to small emitting 
areas. However, another way of opening gaps in disks is 
to expose the disks to energetic radiation. Here we dis- 
cuss the viability of photoevaporation as a hole-opening 
mechanism. We show that, although photoevaporation 
can be responsible for transitional disk holes about 3- 
6 AU wide, it cannot build the extremely wide holes that 
concern us here. 

X-ray, extreme-UV and far-UV photons transfer ki- 
netic energy to gas in the disk, either directly or through 
the photoelectric effect, allowing the gas to become un- 
bound from the star. Photoevaporative flows begin at 
the graviational radius, where the kinetic energy of the 
radiatively accelerated gas particle is equal to the escape 
speed from the star's potential: 



In Equationll4| Cg is the sound speed in the heated upper 
layers of the disk. The localization of the photoevapora- 
tive flow lends itself well to opening gaps in disks, which 
may then widen into holes as the interior disk drains 
on the viscous timescale. However, recall the discussion 
of the viscous timescale in Sj3j in order for accretion to 
drain away interior disk material within the lifetime of 
the disk, the gap must be located within a few AU of the 
star. Of X-ray, EUV and FUV photons, only the X-rays 
have a gravitational radius that might meet our diffusion 
requirements. 

In X-ray photoevaporation, disk clearing begins at a 
characteristic "clearing age" of 2-4 Myr (Owen et al. 
2010). A gap in the disk opens at a characteristic radius 
of 3 AU and the outer and interior disks are decoupled, 
so that the interior disk drains on the viscous timescale, 
several 10^ years to become optically thin for the surface 
density profile in Equation |8l X-rays take approximately 
10-20% of the clearing age, or 2 x 10^-8 x 10^ years, 
to deplete the material between 3 and 10 AU. Viscous 
accretion of material inside 3 AU, to leave behind only 
optically thin dust, and outward expansion of the gap 
therefore occur on roughly the same timescale. Once the 
interior disk has been accreted, the disk outside 10 AU is 
eroded by the X-ray field hitting the hole's inner wall di- 
rectly on a timescale of order the clearing age, a few mil- 
lion years. Accreting transitional disks with inner holes 
between 3 and 10 AU wide are therefore consistent with 
the X-ray photoevaporation model. Indeed, Owen et al. 
(2010) found that approximately 50% of observed transi- 
tion disks could be explained by X-ray photoevaporation. 
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A closer look at the X-ray photoevaporation process, 
however, reveals why X-ray driven winds cannot explain 
our test sample of transitional disks. Once clearing be- 
gins to proceed from 10 AU outward via the direct X-ray 
field, there is no longer any interior mass reservoir to pro- 
vide accretion onto the star. Accreting disks with large 
inner holes therefore cannot be explained by X-ray driven 
winds. For a 20 AU inner hole, such as in the GM Aur 
disk, Owen et al. (2010) find a probability of less than 
2 X 10"'* that the star should still be accreting (see their 
Figure 12). Furthermore, according to the Owen et al. 
models, one would have to wait at least 4 Myr to get a 
20 AU hole — 2 Myr at minimum for the initial gap to 
open and another 2 Myr for it to move out to 20 AU. 
Their models predict that transitional disks with wide 
holes should be older than the average disk population, 
whereas Furlan et al. (2009) and Merin et al. (2010) ob- 
serve no correlation between transitional disk fraction 
and cluster age. 

We next consider photoevaporation driven by EUV ra- 
diation from either the central star or nearby, massive 
stars. EUV continuum radiation photoionizes the upper 
layers of a disk so that hydrogen recombination creates 
a diffuse radiation field that steadily evaporates H II gas 
from the disk surface (HoUenbach et al. 1994, Clarke et 
al. 2001). Once the inner disk has become optically thin 
to Lyman continuum radiation, the direct UV field also 
drives a photoevaporative wind that drains the outer disk 
(Alexander et al. 2006a). However, extreme UV photo- 
evaporation only drives disk dissipation when the vis- 
cous mass flow through the disk becomes weaker than 
the photoevaporative flow, estimated by Gorti & HoUen- 
bach (2009) as 

Msc/y < IO-'-^Mq yr-i. (15) 

The stars hosting transitional disks in our test sample 
have higher accretion rates (Table [T]). Models of UV 
photoevaporation therefore indicate that it cannot be re- 
sponsible for clearing the holes in our test sample. 

Finally, we consider FUV-induced photoevaporation, 
for which the gravitational radius is outside 30 AU (Mat- 
suyama et al. 2003). Such a large gravitational radius 
is appealing as a mechanism of opening wide holes, but 
once more, we have a timescale problem with clearing the 
disk interior to the gravitational radius. For our recon- 
structed GM Aur surface density profile (Eq. [8|) , the sin- 
gle e- folding diffusion timescale from 30 AU is 1.7 Myr- 
higher than the age of GM Aur and other transitional 
disks in Taurus. Furthermore, because the FUV-driven 
wind outside the gravitational radius is so much more ef- 
ficient than viscous accretion inside the gravitational ra- 
dius, FUV-induced photoevaporation tends to truncate 
disks rather than opening inner holes (Matsuyama et al. 
2003). 

Our final piece of evidence that photoevaporation is 
not responsible for the holes in our sample of transitional 
disks (Table[T]) is its inability to create non-axisymmetric 
structure: photoevaporative clearing proceeds evenly 
outward from the gravitational radius. Photoevapora- 
tion therefore does not set up the confined accretion flow 
necessary to move mass efficiently through an optically 
thin hole. Multiple planets provide the only disk-clearing 
mechanism that can simultaneously account for the wide 



inner holes, high accretion rates, confinement of hot gas, 
and asymmetry observed in transitional disks. 

7. CONCLUSIONS AND TESTABLE PREDICTIONS 

While the initial discovery of transitional disks with 
wide holes were based on near- and niid-IR fiux deficits 
(often measured by the 13 /xm/30 fj,m flux ratio; Furlan 
et al. 2006), more detailed observational follow-up of a 
subset of these disks has revealed a number of unusual 
characteristics, including: 

1. Holes of 15 AU~70 AU or more that are too wide 
to have been cleared by a single planet (Brown et 
al. 2007, Brown et al. 2009, Hughes et al. 2010, 
Espaillat et al. 2010, Andrews et al. 2011); 

2. Stellar- mass companions orbiting within 4 AU of 
the central star ruled out at the 99.7% confidence 
level for three of our disks in Table [T] (Pott et al. 
2010); 

3. Accretion onto the star at rates 10^^ < M < 
10-8 Mq yr-i (White & Ghez 2001, Johns-Kruh 
& Gafford 2002, Andrews et al. 2011); 

4. Complex radial mass distributions of dust and 
gas derived from SED models, millimeter im ages, 
IR veiling, and IR interferometry jE spaillat etaLl 
[2OIOI: lAndrews et al.|[20Tlt lAkeson et al...20ll : 

5. Non-axisymmetric structure derived from millime- 
ter images , spectro-astr o metry and sparse aperture 
mask ing ("Brown et a l.l 120091 iPontoppidan et aD 
120081 : E^elamo ct al. 20 lit). 

In this article, we have shown that multiple plane- 
tary systems are the only phenomena that can explain 
all of the transitional disk characteristics listed above. 
The other possible disk-clearing mechanisms — massive 
single planet, companion star and photoevaporation — 
have timescale problems for wide holes and do not pro- 
vide ways to hide accreting mass in holes that appear 
optically thin. 

The planetary system configuration in the multiplanet 
model must satisfy three constraints. First, each planet 
must be a gas giant of ^ IMj or more so that it can open 
a gap. Second, the planets must remain on stable orbits 
for the lifetime and/or clearing time of the transitional 
disk lest they disrupt the outer, optically thick disk out- 
side the hole. Finally, for sources with extremely low 
near-IR fluxes such as GM Aur, the innermost planet 
must be located within ~ 2 AU of the central star so 
that the interior disk pileup can diffuse away within the 
hole-clearing timescale. For the disks such as DoAr 44 
with large holes in millimeter images but higher near-IR 
fluxes, the constraint on the location of the innermost 
planet is somewhat relaxed. 

Central to our success at finding a single, unifying the- 
ory for the wide-hole transitional disks is our break away 
from assuming axisymmetric, radially uniform dust and 
gas distributions. In illustrating how accreting mass can 
be hidden in apparently optically thin holes (Figures 1, 
[3] and [5]) , we are using a geometrical definition of optical 
depth — the macroscopic fraction of an emitting surface 
area that is covered with light sources. Tidal tails and 
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narrow annuli of dust and gas trapped between planets 
are locally optically thick, but cover only a small frac- 
tion of the inner hole. Overall, then, the inner hole can 
appear optically thin, but still transport mass onto the 
star. Our hydrodynamic simulations suggest that taking 
into account the possibility of a non-uniform dust and 
gas distribution at least doubles the amount of dust that 
is consistent with an optically thin inner hole. 

The assumption of axisymmetrically distributed gas 
and dust is problematic because it may hide the dy- 
namical effects of planetary systems. The existence of 
confined "funnel flows" raises the possibility that time 
variability of the star accretion rate or excess near-IR 
flux may be a signature of planetary influence. Accre- 
tion rate or near-IR variability would occ ur on the orbita l 
timescale of the innermost planet. Espail lat et al.l ()2011l ) 
detected variability on 2-3 year timescales in all but two 
of the 14 transitional and pre-transitional disks in their 
study of the Taurus and Chamaeleon star-forming re- 
gions. 

The level of axisymmetry in inner disk flow patterns is 
an observational test of our multiplanet theory of transi- 
tional disk origins. As discussed in fJJJ non-axisymmetric 
structures h ave been discovered in a number of transi- 
tional disks (iBrown et al.ll2008l: Fontoppidan et al.ll2008l: 
iHuelamo et al.l 1201 1[ ). Embedded protoplanets are also 
expected t o produce asy mmetries in emission line pro- 
files (Rcg alvet al.l [20Tol ). Furthermore, the finding of 
Salyk et al. (2011) that hot CO gas in transitional disks 
has a different, possibly more confined, radial profile as 
compared to classical disks should be investigated fur- 
ther with disk models. If confirmed, hot gas confinement 
strongly indicates a dynamical origin for transitional disk 
holes rather than radiation-driven clearing. 

Finally, we offer a possible explanation, based on 
our multiplanet theory, for the intriguing SR 21 disk 
which we mentioned in ^ SR 21 has a gas disk 
with a 7 AU inner hole but a dust disk with a 
30 AU hole fPonto ppidan et al.ll2008l: IBrown et al.]|200l 
lAndrews e t al. 20111). If grain growth and planetesimal 
formation in SR 21 have proceeded to the point where 
very few primordial, small grains remain, any observed 
dust would be second-generation, supplied by colliding 
planetesimals. As we have seen in ij3.1l giant planets are 
extremely effective at scattering test particles away from 
their orbits. A system of gas giants extending to ^ 30 AU 
would deplete the region of planetesimals, choking off any 
chance of producing second-generation dust. Gas, how- 
ever, could still accrete through the planetary system, 
explaining the smaller gas truncation radius. The SR 21 
system fits with our multiplanet theory as long as the 
observed dust is primarily the product of planetesimal 
collisions. 

As high angular-resolution astronomy advances, we 
urge observers to focus planet-hunting campaigns on 



transitional disks. ALMA will provide detailed kinematic 
and spatial information about the gas and dust distribu- 
tions in transitional disk holes. For the extremely wide, 
70 AU holes reported by Andrews et al. (2011), direct 
imaging with GPI may yield pictures of young planets 
in formation. Finally, aperture-masking techniques have 
the potential to reveal planetary-m ass objects within 
15 AU of their central stars (e.g. iKraus et al.ll2009l) 
and may uncover the hidden planetary systems we sug- 
gest are inside the transitional disk holes. In fact, such a 
companion may have already been detected in the wide 
transitional disk T Cha (jHuelamo et al.ll20lil) . Intrigu- 
ingly, the only planets that can open gaps in disks are 
gas giants of order Jupiter's mass. Transitional disks 
have excellent discovery potential for young, still-forming 
planets, as planet searches focusing on young stars are 
sensitive only to massive planets. 

The discovery of extremely wide holes in transitional 
disks (Brown et al. 2007, Andrews et al. 2011) points to 
a possible gravitational-instability origin for the planets 
contained in the holes: Dodson-Robinson et al. (2009b) 
demonstrated that gas giants with semimajor axes > 
35 AU can only form by gravitational instability. As 
a bonus, gravitational instability tends to produce hot- 
ter, more massive planets than core accretion (Clarke 
2009, Stamatellos & Whitworth 2009), exactly the types 
of planets best suited to discovery by direct imaging and 
aperture masking. In contrast, radial-velocity observa- 
tions of main-sequence stars, which have provided almost 
all of the 500 and more exoplanets discovered to date, are 
strongly biased toward planets formed by core accretion, 
which is effective for shorter-period orbits. Short-period 
planets' inability to open wide holes in disks could also 
explain why radial velocity surveys of young stars with 
disks, which were des igned to include l arge fractions of 
transitional disks (e.g. lPrato et al.l[2008( ). have yet to de- 
tect a hot Jupiter. In any case, studies of transitional 
disks will likely provide insight into a newly confirmed 
and only partially understood planet formation mecha- 
nism. 
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APPENDIX 



PROTOSTELLAR DISK DIFFUSION TIMESCALE 



To calculate the rate at which disk material accretes onto the star, we consider the surface density diffusion equation 
(Pringle 1981): 



5S 
dt 



3d_ 

r dr 



,.1/2 A 

dr 
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in which S is the surface density of disk material, r is the distance from the star, t is time, and v is the kinematic 
viscosity of the disk material Over a small region of the disk, one can approximate the diffusion rate by assuming a 
constant viscosity so that the diffusion equation can be separated into a time component and a space component: 

nr.t)=R{r)T{t). (A2) 

iMassetl (|200Clf) and iMasset fc Papaloizoul ()2003[ ) also used disks with constant kinematic viscosity in their analytical 
and numerical studies of planet migration. By taking the partial derivatives of the separated function I](r, i), we 
rewrite the diffusion equation as 

1 T' 9 i?' i?" , , , , 

where T' and R' are partial derivatives over t and r, respectively, and A is the constant Eigenvalue. The solution to 
the time-dcpcndcnt portion of equation IA3I is 

T{t) = exp {-\vt) . (A4) 

To find an expression for A that is benchmarked against measured T-Tauri star accretion rates and does not depend 
on assumptions about viscosity, we examine the mass flow through the interior disk. For an interior disk annulus of 
finite size such as the one left behind by a gap-clearing planet, one may approximate the total mass transfer through 
the annulus, M, as 

-M^i:rlt{rg), (A5) 

where Tg is the inner edge of the gap created by the planet. The negative sign reflects the convention that material 
flowing inward creates a positive accretion rate onto the star. Substituting E(r) = R{r)T'{t) into equation IA5I and 
solving for A, we find 

M 

A«^^, (A6) 

where E is the average surface density in the diffusing region r < Vg. 
Using Equation I A6 1 we may write the solution to the diffusion equation under the assumption of constant viscosity 

as 

E(r,t)«S(r,0)exp|^^^y (A7) 

where Mq is the accretion rate at t = 0. Here the diffusion timescale we recover is similar to the canonical value 
^diff = r^/i/ (Lynden-Bell & Pringle 1974), but in terms of the accretion rate: 

idiff = (A8) 

Mq 

The quantity TrrgE is, of course, equivalent to the disk mass Mivg) contained inside a circle of radius r^, so we rewrite 
Equation I A8I to get Equation [7] in fJH 

For disks with E cx r~^, as is usually inferred from sub-mm observations, the diffusion timescale from Equation [7] 
is linearly proportional to r^. A standard a- viscosity disk model with steady-state mass transfer has u oc r so that 
the radial dependence of idiff according to the traditional /v criterion is also linear. Even though we have, for the 

purpose of clarity, based our derivation of Equation [7] on a disk with constant viscosity, the linear dependence of M/M* 
on radius applies to a range of commonly used disk models. 

ACCRETION EFFICIENCY AND PLANETARY INTERIOR STRUCTURE 

Bodenheimer (1974) divides the evolution of a young gas giant into four phases. The first is a transient hydrodynamic 
epoch, in which gas at the planet's Hill radius free-falls toward the center of the planet: the planet's rapid growth is 
driven by the runaway expansion of its Hill sphere. The hydrodynamic epoch lasts about one free-fall time, ~ 1,000 
years. For planets growing by core accretion, Phase I encompasses the mass range of ^ lO-SOOM^. For planets formed 
by disk fragmentation, the characteristic Phase I mass is much higher (e.g. Clarke 2009). Planets in Phase I accrete 
100% of the material inside their Hill spheres. Gap opening begins during Phase I at about Saturn's mass, 95M0, and 
eventually halts the runaway growth of the planet, which is why we do not start our simulations with small planet 
cores and let them accrete: FARGO 's accretion efficiency parameter p would suffer a dramatic drop during only the 
first 1,000 years of disk evolution. 

At the onset of Phase H, the young planet develops an interior pressure gradient that halts the free-fall of material at 
the Hill radius. The planet then stays in Phase H, quasistatic contraction on the Kelvin-Helmholz timescale over about 
a factor of 20 in radius, for ^ 10^ years. The planet's evolution during this period is similar to that of a protostar on 
the Hayashi track. Accretion during Phase H is limited by the disk's ability to deliver gas to the edge of the gap and 
by the strength of the planet's tidal barrier (e.g. Dobbs-Dixon et al. 2008). "Leakage" of gas across the tidal barrier 
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is much more cfBcicnt at the beginning of Phase II, when the planet nearly fills its Hill sphere: at this stage, material 
from the disk can fall directly to the planet's surface. As the planet contracts, the velocity of disk material entering the 
Hill sphere increases relative to the planetary surface, so that accreting gas will orbit the planet in a circumplanetary 
disk (Papaloizou & Nelson 2005). Mass transfer onto the planet is then regulated by the the diffusion timescale in the 
circumplanetary disk as well as the mass transfer rate through the protostellar disk and is extremely inefficient (Ward 
& Canup 2010). 

Once the planet's interior reaches 2000 K, molecular hydrogen begins to dissociate, draining the thermal energy that 
maintains the quasistatic pressure gradient. The resulting supersonic collapse, which starts Phase HI, takes only a few 
days and decreases the planet's radius by a factor of more than 100, further limiting the planet's ability to accrete. 
We see that a young planet's accretion efficiency changes dramatically as it evolves. During most of the transitional 
disk epoch, the planet will be accreting through a circumplanetary disk whose mass transfer efficiency depends on the 
ratio of protostellar/circumplanetary disk viscosities and the size of the planet. The circumplanetary disk viscosity 
depends in turn on the density and temperatme of circumplanetary material. 

Since FARGO is an isothermal code that docis not treat the planetary interior structure, it much better suited for 
parameter studies of gap-opening efficiency in the presence of planetary accretion than physically realistic computations 
of planet growth profiles. Because the accretion-limiting circumplanetary disk should be present and growing during 
most of the hole-clearing epoch, we are justified in assuming that mass transfer onto the planets is negligible throughout 
our FARGO simulation. 



14 



Dodson-Robinson & Salyk 



REFERENCES 



Akeson, R. L., et al. 2005, ApJ, 635,173 
Akeson, R. L., et al. 2011, ApJ, 728, 96 

Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 
369, 216 

Alexander, R. D., & Armitage, P. J. 2007, MNRAS, 375, 500 
Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C, & 

DuUemond, C. P. 2010, ApJ, 723, 1241 
Andrews, S., Wilner, D., Espaillat, C, Hughes, A. M., 

DuUemond, C. P., McClure, M. K., Qi, C, & Brown, J. M. 

2011, ApJ, in press 
Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 
Balbus, S. A., Hawley, J. F., & Stone, J. M. 1996, ApJ, 467, 76 
Bast, J. E., Brown, J. M., Herczeg, G. J., van Dishoeck, E. F., & 

Pontoppidan, K. M. 2011, A&A, 571, 119 
Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. 

1990, AJ, 99, 924 
Bodenheimer, P. 1974, Icarus, 23, 319 

Brown, J. M., Blake, G. A., DuUemond, C. P., Merfn, B., 
Augereau, J. C., Boogert, A. C. A., Evans, N. J., H, Geers, 
V. C, Lahuis, F., Kessler-Silacci, J. E., Pontoppidan, K. M., & 
van Dishoeck, E. F. 2007, ApJ, 664, L107 

Brown, J. M., Blake, G. A., Qi, C, DuUemond, C. P., & Wilner, 
D. J. 2008, ApJ, 675, L109 

Brown, J. M., Blake, G. A., Qi, C, DuUemond, C. P., Wilner, 
D. J., & Williams, J. P. 2009, ApJ, 704, 

Calvet, N., D'Alessio, P., Hartmann, L., Wilner, D., Walsh, A., & 
Sitko, M. 2002, ApJ, 568, 1008 

Calvet, N., et al. 2005, ApJ, 630, L185 

Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, 
ApJ, 686, 580 

Chiang, E., & Murray-Clay, R. 2007, Nature Physics, 3, 604 
Clarke, C. J. 2009, MNRAS, 396, 1066 

Clarke, C. J., Gendrin, A., & Sotomayor, M. 2001, MNRAS, 328, 
485 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 
245 

Crida, A., Morbidelli, A., & Masset, F. 2006, Icarus, 181, 587 
Currie, T., & Sicilia-Aguilar, A. 2011, ApJ, in press 

(arXiv: 1102.4264) 
Dobbs-Dixon, I., Li, S., & Lin, D. N. C. 2008, lAU Symposium, 

249, 263 

Dodson-Robinson, S. E., Veras, D., Ford, E. B., & Beichman, 

C. A. 2009, ApJ, 707, 79 
Dodson-Robinson, S. E., Willacy, K., Bodenheimer, P., Turner, 

N. J., & Beichman, C. A. 2009, Icarus, 200, 672 
Draine, B. T., et al. 2007, ApJ, 663, 866 
DuUemond, C. P., & Dominik, C. 2004, A&A, 417, 159 
Edwards, S., Fischer, W., Hillenbrand, L., & Kwan, J. 2006, ApJ, 

646, 319 

Eisner, J. A, Chiang, E. I., & Hillenbrand, L. A. 2006, ApJ, 637, 
L133 

Ercolano, B., Drake, J. J., Raymond, J. C, & Clarke, C. C. 2008, 
ApJ, 688, 398 

Espaillat, C, Calvet, N., D'Alessio, P., Hernandez, J., Qi, C, 
Hartmann, L., Furlan, E., & Watson, D. M. 2007, ApJ, 671, 
L135 

Espaillat, C, Calvet, N., Luhman, K. L., MuzeroUe, J., & 

D'Alessio, P. 2008, ApJ, 682, L125 
Espaillat, C, et al. 2010, ApJ, 717, 441 

Espaillat, C, Furlan, E., D'Alessio, P., Sargent, B., Nagel, E., 
Calvet, N., Watson, D. M., & MuzeroUe, J. 2011, ApJ, 728, 49 

Fernandez, M., Ortiz, E., Eiroa, C, & Miranda, L. F. 1995, 
A&AS, 114, 439 

Folha, D. F. M., & Emerson, J. P. 1999, A&A, 352, 517 

Furlan, E., Hartmann, L., Calvet, N., D'Alessio, P., 

Franco-Hernandez, R., Forrest, W. J., Watson, D. M., Uchida, 
K. I., Sargent, B., Green, J. D., Keller, L. D., & Herter, T. L. 
2006, ApJS, 165, 568 

Furlan, E., Watson, D. M., McClure, M. K., Manoj, P., Espaillat, 
C, D'Alessio, P., Calvet, N., Kim, K. H., Sargent, B. A., 
Forrest, W. J., & Hartmann, L. 2009, ApJ, 703, 1964 

Garcia Lopez, R., Natta, A., Testi, L., & Habart, E. 2006, A&A, 
459, 837 

Goldreich, P., & Tremaine, S. 1982, Annual Review of Astronomy 
and Astrophysics, 20, 249 



Gorti, U., & HoUenbach, D. 2009, ApJ, 690, 1539 

Hartmann, L., Megeath, S. T., Allen, L., Luhman, K., Calvet, N., 

D'Alessio, P., Franco-Hernandez, R., & Fazio, G. 2005, ApJ, 

629, 88 

Hatzes, A. P., et al. 2000, ApJ, 544, L145 

Hawley, J. F., & Stone, J. M. 1998, ApJ, 501, 758 

Herbst, W., Herbst, D. K., Grossman, E. J., & Weinstein, D. 

1994, AJ, 108, 1906 
HoUenbach, D., Johnstone, D., Lizano, S., & Shu, F. 1994, ApJ, 

428, 654 

Huelamo, N., Lacour, S., Tuthill, P., Ireland, M., Kraus, A., & 

Chauvin, G. 2011, A&A, 528, L7 
Hughes, A. M., Andrews, S. M., Espaillat, C, Wilner, D. J., 

Calvet, N., D'Alessio, P., Qi, C, Williams, J. P., & 

Hogerheijde, M. R. 2009, ApJ, 698, 131 
Ireland, M. J., & Kraus, A. L. 2008, ApJ, 678, L59 
Johns-KruU, C. M., & Gafford, A. D. 2002, ApJ, 573, 685 
Kalas, P., Graham, J. R., Chiang, E., Fitzgerald, M. P., Clampin, 

M., Kite, E. S., Stapelfeldt, K., Marois, C, & Krist, J. 2008, 

Science, 322, 1345 
Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 
Kley, W., Bitsch, B., & Klahr, H. 2009, A&A, 506, 971 
Kley, W., & Crida, A. 2007, A&A, 487, 9 
Koerner, D. W., Sargent, A. I., & Beckwith, S. V. W. 1993, 

Icarus, 106, 2 

Kraus, A. L., Ireland, M. J., Martinache, F., & Lloyd, J. P. 2009, 

American Institute of Physics Conference Series, 1094, 453 
Kraus, A. L., Ireland, M. J., Martinache, F., & Hillenbrand, L. A. 

2011, ApJ, 731, 8 
Lagrange, A.-M., Bonnefoy, M., Chauvin, G., Apai, D., 

Ehrenreich, D., Boccaletti, A., Gratadour, D., Rouan, D., 

Mouillet, D., Lacour, S., & Kasper, M. 2010, Science, 329, L57 
Li, A., & Draine, B. T. 2001, ApJ, 554, 778 
Lin, D. N. C, & Papaloizou, J. C. B. 1993. Protostars and 

Planets III, 749, 835 
Lynden-Bell, D., & Pringle, J. E. 1974, MNRAS, 168, 603 
Marois, C, Macintosh, B., Barman, T., Zuckerman, B., Song, I., 

Patience, J., Lafreniere, D., & Doyon, R. 2008, Science, 322, 

1348 

Marois, C, Zuckerman, B., Konopacky, Q. M., Macintosh, B., & 

Barman, T. 2010, Nature, 468, 1080 
Masset, F. 2000, A&AS, 141, 165 

Masset, F. S., & Papaloizou, J. C. B. 2003, ApJ, 588, 494 
Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 
Matsuyama, I., Johnstone, D., & Hartmann, L. 2003, ApJ, 582, 
893 

Merfn, B., et al. 2010, ApJ, 718, 1200 

MuzeroUe, J., Allen, L. E., Megeath, S. T., Hernandez, J., & 

Gutermuth, R. A. 2010, ApJ, 708, 1107 
Natta, A., Testi, L., MuzeroUe, J., Randich, S., Comeron, F., & 

Persi, P. 2004, A&A, 424, 603 
Natta, A., Testi, L., & Randich, S. 2006, A&A, 452, 245 
Ossenkopf, V., Henning, T., & Mathis, J. S. 1992, A&A, 261, 567 
Owen, J. E., Ercolano, B., & Clarke, C. J. 2011, MNRAS, 210 
Paardekooper, S.-J., Baruteau, C, Crida, A., & Kley, W. 2010, 

MNRAS, 401, 1950 
Paardekooper, S.-J., Baruteau, C, & Kley W. 2011, MNRAS, 

410, 293 

Paardekooper, S.-J., & Mellema, G. 2006, A&A, 450, 1203 
Palla, F., & Stabler, S. W. 2000, ApJ, 540, 255 
Papaloizou, J. C. B., & Nelson, R. P. 2005, A&A, 433, 247 
Pontoppidan, K. M., Blake, G. A., van Dishoeck, E. F., Smette, 

A., Ireland, M. J., & Brown, J. 2008, ApJ, 684, 1323 
Pontoppidan, K. M., Salyk, C, Blake, G. A., Meijerink, R., Carr, 

J. S., & Najita, J. 2010, ApJ, 720, 887 
Pott, J.-U., Perrin, M. D., Furlan, E., Ghez, A. M., Herbst, 

T. M., & Metchev, S. 2010, ApJ, 710, 265 
Prato, L., Huerta, M., Johns-KruU, C. M., Mahmud, N., Jaffe, 

D. T., & Hartigan, P. 2008, ApJ, 687, L103 
Pringle, J. E. 1981, ARA&A, 19, 137 

Quillen, A. C, Blackman, E. G., Frank, A., & Varniere, P. 2004, 

ApJ, 612, L137 
Regaly, Zs., Sandor, Zs., DuUemond, C. P., & van Boekel, R. 

2010, A&A, 523, 69 



Young multiplanet systems 



15 



Rcttig, T. W., Haywood, J., Simon, T., Brittain, S. D., & Gibb, 

E. 2004, ApJ, 616, L163 
Rice, W. K. M., Armitage, P. J., Wood, K., & Lodato, G. 2006, 

MNRAS, 373, 1619 
Rice, W. K. M., Wood, K., Armitage, R J., Whitney, B. A., & 

Bjorkman, J. E. 2003, MNRAS, 342, 79 
Robbins, S. J., Stewart, G. R., Lewis, Mark C., Colwell, J. E., & 

Sremcevic, M. 2010, Icarus, 206, 431 
Salyk, C., Blake, G. A., Boogert, A. C. A., & Brown, J. M. 2009, 

ApJ, 699, 330 
Salyk, C., et al. 2011, in prep 

Schisano, E., Covino, E., Alcala, J. M., Esposito, M., Gandolfi, 

D., & Guenther, E. W. 2009, A&A, 501, 1013 
Scharf, C., & Menou, K. 2009, ApJ, 693, 113 
Stamatellos, D., & Whitworth, A. P. 2009, MNRAS, 392, 413 
Tanaka, H., Takeuchi, T., & Ward, W. 2002, ApJ, 565, 1257 
Turner, N. J., & Drake, J. F. 2009, ApJ, 703, 2152 



Turner, N. J., & Sano, T. 2008, ApJ, 679, 131 
Varnicrc, P., Blackman, E. G., Prank, A., & Quillen, A. C. 2006, 
ApJ, 640, 1110 

Ward, W. R. 1991, Lunar and Planetary Science Institute Science 

Conference Abstracts, 22, 1463 
Ward, W. R. 1997, Icarus, 126, 261 
Ward, W. R., & Canup, R. M. 2010, AJ, 140, 1168 
Weaver, W. B., & Jones, G. 1992, ApJS, 78, 239 
Weidenschilling, S. J. 1977, MNRAS, 180, 57 
Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296 
White, R. J., & Ghez, A. M. 2001, ApJ, 556, 265 
Zhu, Z., Nelson, R. P., Hartmann, L., Espaillat, C., & Calvet, N. 

2011, ApJ, 729, 47 
Zubko, V. G., Mennella, V., ColangeU, L., & Bussoletti, E. 1996, 

MNRAS, 282, 1321 



16 



Dodson-Robinson & Salyk 




Wavelength [jum] Radius [AU] 

Fig. 1. — In this figure we demonstrate some of the degeneracies inherent to modeling of near- and mid-IR excess flux, especially once 
radial and azimuthal asymmetry are permitted. The left panel shows near- and mid-IR excess fluxes measured for GM Aur with two model 
SEDs based on different mass distributions (square: Edwards et al. 2006; open circles: Foiha & Emerson 1999 Q ; filled circles: Espaillat et 
al. 2010 E; star: Wea ver & Jones 1992; spectru m: Calvet et a l. 2005, Pontoppidan et al. 2010). Q Fluxes are de-reddening using Ay = 1.2 
IIEspaillat et al.ll20ig ) and the reddening law of lCardelli et all 11989. ). The model mass distributions are shown schematically, and to scale, 
in the right hand panels, where gray represents optically thin regions, black represents optically thick regions, and white represents empty 
regions. The outer disk contribution from R>20 AU is not shown, but primarily contributes at A > 13/im. 

''These veiling values may be high due to a mismatch in photospheric template, or variability Ij EspaiUat et al."2010'). 

°Last three points are estimated from Figure 5 in lEspaillat eFaLl | |20TCD and are assigned the er ror bar a ssociated with rj^. 

'^We omit the M-band veiling values because of a likely template mismatch in their calculation IjSalvk et al.ll2009f l. 
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Orbital Radius (AU) 

Fic;. 2. — The scattering timcscale of particles initially on circular orbits approaches the protostellar disk lifetime for particles more than 
~ 5 Hill radii from a planet. Here we plot scattering tiniescale as a function of distance from the star in a disk containing a planet at 10 AU 
(vertical dotted line). The high scattering timescales interior to the planet's orbit demonstrate that some other disk dissipation mechanism 
must cLct along with planetary clearing to clean out the inner holes of transitional disks. 
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Fig. 3. — Top Row; At 20 AU wide, the hole in the GM Aur disk is too large to have been cleared by a single planet. Three medium- 
mass planets, all 3A/j with distances 14.3 AU, 6.3 AU and 2.7 AU from GM Aur, clear gaps that overlap and form a large hole over the 
course of 300 orbits of the outer planet (2 X 10* years). Here we show plots of the disk surface density from orbits 2, 197, and 395 of the 
FARGO simulation, with the color scale chosen so white regions are optically thin at millimeter wavelengths. The axes show position in the 
corotating frame, where the outermost planet is at (x,y) = (1,0). Bottom Row; Here we demonstrate that tidal streams can efficiently 
transport mass through a hole that appears to be optically thin in the near- and mid-infrared. Simulation timepoints are the same as in 
the top row, but the color scale is designed so that all white regions are optically thin at 13 pixa — assuming gas/solid mass ratio of 100, 
all of the disk's solid mass is in micron or submicron grains, and k;^ dust = 500 cm-^ g~^ for the grains {Osscnkopf et al. 1992). Note how 
the non-axisymmetric, confined flow pattern produced by the planets keeps most of the hole empty while allowing mass transport through 
locally optically thick tidal streams. 
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Fig. 4. — Top Row: Renderings of the 400 fim optical depth of our model disk at orbits 197 and 395, here assuming that 80% of the 
disk's solid mass is in planetesimals and planet cores that do not contribute to millimeter flux. As in Figure|3l all white regions are optically 
thin. Bottom Row: The 13 fim optical depth of our model disk if 80% of the solids have formed pebbles, planetesimals or planets. After 
400 orbits, locally optically thick dust is confined to only about 10% of the hole surface area, yet the inner hole contains about Saturn's 
mass in gas and solids. Here we assume an opacity of kx = 500 cm-^g~^ for micron/submicron grains. The total gas/solid mass ratio is 
still 100. 
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Fig. 5. — Left: While the azimuthally-averaged dust surface density S^j in our model transitional disk (black solid line) implies optically 
thick dust everywhere except in two narrow annuli near the outermost planet, taking into account azimuthal asymmetries yields a filling 
factor, /, (gray solid line) less than unity over a substantial portion of the inner disk. The dotted line denotes the optically thin dust surface 
density cutoff at 13 /im. The filling factor f{r) w as co mputed assuming a gas/small grain mass ratio of 100:1, but realistic gas/grain ratios 
in plane t-forming disks may be mu ch higher (see ^4.2l and Figure|4]l. Center: A comparison between the observed SE P of GM Aur, shown 
in gray IIKenvon Hartmann|[T995i : iWeaver fc Joneail992 : iHartmann et al.i,2005i : iCalvet et al.ii2005l : IPontoppidan et al.i 2010). and fiuxes 
produced using our reference multi-planet model. Models consist of an optically thin inner disk (dotted line, but below plot range for GM 
Aur), an optically thick inner disk meant to represent our dynamical model (dashed line shows best-fit case), and an optically thick outer 
disk (dot-dashed line). From top to bottom, models have / = 1 (red squares), / = 0.5 (blue triangles) and / = 0.03 (black). Right: A 
comparison between the observed SED of DoAr 44 (references shown above plus Herbst et al. 1994) and our reference model. From top to 
bottom, models have / = 1 (red squares), / = 0.5 (blue triangles) and / = 0.1 (black). 
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TABLE 1 

Observed hole sizes and accretion rates of 
transitional disks 

Name hole size Ref M Ref 
[log(Ms yr-i)] 



DM Tau 


19 


m 


-8.7 


DoAr 44 


30 


Ihlhl 


-8.0 


GM Aur 


20 


-8.0 


HD 135344 


39 


m 


-8.3 


Lk Ha 330 


68.47 


mm 


-8.8 


MWC 758 


73 




-8.0 


SR 24 S 


29 




-7.2 


T Cha 


15 


m 


-8.4 


UX Tau A 


25 


m 


-8.0 



m 



m^ 



m 



References. — (1) Andrews et al. 2011 (2) Brown 
et al. 2007 (3) Brown et al. 2009 (4) Espaillat et al. 2010 
(5) Fernandez et al. 1995 (6) Garcia Lopez et al. 2006 (7) 
Hughes et al. 2009 (8) Johns-Krull & Gafford 2002 (9) 
Natta et al. 2006 (10) Schisano et al. 2009 (11) White & 
Ghez 2001 

Estimate d from full w idth of Ha at 10% level and Equa- 
tion (1) of lNatta et"all f2004) . 
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TABLE 2 
Parameters of FARGO simulation 



Parameter 



Description 



Value 



Disk setup 

Tin inner boundary 

''out outer boundary 

grid zones 

boundary conditions 
Hjr aspect ratio 
a viscous efficiency 

planet release time'' 
Planet setup 

Mp planet masses 
r-p initial orbit radii*^ 

initial planet separation 



0.72 AU 

38 AU 

220 radial, 384 azimuthal 
open with outer mass source'' 
0.05 
0.002 

12,000 years 
3Mj 

2.7 AU, 6.3 AU, 14.3 AU 



Open boundaries allow mass to move freely through the grid and accrete 
onto the star. The outer mass source allows material from beyond 38 AU, 
which we do not model, to move into the grid assuming steady-state 
accretion. 

To remove unphysical Type I migration torques from our simulation, 
we hold the planets on fixed circular orbits for 12,000 years (200 orbits of 
the outermost planet), then turn on migration and fewbody interactions. 

AU planets were initially placed on circular orbits. 



